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NOMENCLATURE 

b  boundary  layer  parameter 

Cf  friction  coefficient 

c  absolute  velocity  component,  chord 

E  entrainment  rate 

F  force  defect  thickness 

f  blade  force,  wake  model  (eq  (7.4)) 

g  wake  model  (eq  (7.7)) 

H  shape  factor 

* 

H  Head's  shape  factor 

h  annulus  height,  metric  coefficient 

m  meridional  coordinate 

n  coordinate  normal  to  the  end-walls 

p  static  pressure 

R^  radii  of  curvature  (appendix  A) 

r  radius,  recovery  factor 

3  streamwi3e  coordinate,  pitch 

T  temperature 

T  total  temperature 

recovery  temperature 

T*  Eckert  reference  temperature 

t  transverse  coordinate 

t  tip  clearance 

c 

U  wheel  speed 

u  pitchwise  coordinate 

W  velocity  in  the  blade  reference  system 

w  relative  velocity  component 


-  3  - 


x  orthogonal  coordinate 


a  radial  flow  angle 

8  passage  flow  angle 

Y  stagger  angle 

6  physical  boundary  layer  thickness 

* 

6  displacement  thickness 
e  skewing  angle,  diffusion  coefficient 
wall  skewing  angle 
n  u/6 

0  momentum  thickness,  angular  coordinate 

A  correction  momentum  thickness 

v  viscosoty,  eddy  viscosity 
p  density 

o  solidity,  wall  inclination  angle  (fig  3.1) 
t  shear  stress 

to  vorticity  vector 

Subscripts 


k 

kinematic 

w 

at  the  wall 

m,n  ,u 

in  meridional  coordinates 

(fig.  3-D 

s,n,t 

in  streamline  coordinates 

(fig.  6.1) 

sec 

secondary  flow 

ewbl 

end-wall  boundary  layer 

inv 

inviscid 

Superscripts 


inviscid  flow  parameter 

passage-averaged  value 

expressed  in  s,t  coordinates  (figure  6.1) 

Abbreviations 


EWBL  end-wall  boundary  layer 
PBL  profile  boundary  layer 
3D  three-dimensional 


1 .  INTRODUCTION 


The  main  objective  of  the  present  research  is  to  provide  a  method  for  the 
prediction  of  detailed  flow  properties  in  a  multistage  axial  compressor,  in 
particular  the  secondary  flow  patterns  and  the  radial  mixing  process. 

The  background  of  the  present  work  lies  in  a  Quasi  3D  approach  [1].  In 
this  approach,  correlations  are  used  to  compute  the  passage-averaged  profile 
losses  and  deviations.  End-wall  boundary  layers  (EWBL)  are  predicted 
through  an  3D  integral  method  where  end-wall  secondary  flows  and  clearance 
effects  are  introduced  [2,3]. 

In  this  EWBL  approach,  it  is  assumed  that  all  the  secondary  and  tip 
clearance  effects  originate  inside  the  EWBL  and  extra  losses  are  taken  into 
account  separately  when  computing  the  machine  efficiency  [3].  However,  the 
interaction  between  the  secondary  flows  and  the  main  compressor  flow  is 
stronger  than  a  simple  efficiency  correction  [4, 5, 6, 7, 8].  The  present  trend 
towards  higher  stage  loadings  and  lower  aspect  ratio's  enhances  the  effect 
of  radial  redistributions  but  current  axisymmetric  calculation  methods  are 
not  able  to  account  for  it.  This  interaction  is  illustrated  on  figure  1.1 
where  a  radial  distribution  of  temperature  is  shown  with  and  without  radial 
redistribution  effects. 

In  the  present  project,  a  technique  for  the  determination  of  the  radial 
mixing  is  being  developed.  This  technique  is  based  on  the  determination  of 
the  radial  flow  components  as  obtained  from  integral  boundary  layer  and  wake 
predictions. 
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Figure  1 . 1  :  Radial  red istr ibut ion  of  total  temperature  [3] 

Adkins  &  Smith  [8]  presented  a  basic  formulation  for  the  determination  of 
radial  mixing  effects.  In  their  approach,  a  large  number  of  parameters 
determining  the  secondary  flow  are  to  be  chosen  on  an  empirical  basis.  In 
the  present  approach,  the  velocity  profile  families  built  in  the  integral 
boundary  layers  and  wakes  lead  to  a  prediction  of  most  of  these  parameters 
and  it  allows  the  reconstruction  of  the  secondary  flow  pattern,  with  muon 
less  additional  parameters. 

As  shown  by  Adkins  and  Smith  [8],  the  radial  mixing  can  be  express^: 
through  a  diffusion  type  equation  where  the  'diffusion'  or  'mixing' 


coefficient  is  a  function  of  the  rms  radial  velocity.  When  the  radial 
velocity  is  well  known  this  mixing  coefficient  can  be  easily  estimated. 
Gallimore  and  Cumpsty  [6]  showed  that  mixing  due  to  turbulent  diffusion  may 
be  as  important  as  mixing  due  to  radial  convection  and  an  additional 
turbulent  mixing  coefficient  is  to  be  considered.  'When  the  total  mixing 
coefficient  is  known,  the  mixing  equation  can  easily  be  solved  using  a 
finite  difference  or  finite  element  method  with  Neumann  type  boundary 
conditions. 


2.  HISTORICAL  BACKGROUND 


During  the  last  twenty  years,  the  design  requirements  for  axial 
compressors  have  increased  considerably  calling  for  more  advanced 
computational  methods.  However,  because  of  the  high  complexity  of  the  flow 
pattern,  no  complete  three-dimensional  calculation  methods  can  presently  be 
developed  which  would  account  for  multistage  configurations.  The  flow  is  to 
be  computed  by  parts  which  interact  with  each  other. 

The  current  approach  consists  in  the  splitting  of  the  turbomachine  flow 
into  two  two-dimensional  flows  [1].  The  first  flow  is  considered  in  a 
circumferential  blade-to-blade  surface  (SI)  and  the  second  is  the 
circumferentially  averaged  meridional  flow  (S2).  Different  secondary 
effects  which  are  not  considered  in  these  SI  and  S2  flows  are  to  be  computed 
separately  and  superposed  in  an  interactive  way.  The  different  secondary 
flows  are  summarised  on  figure  2.1.  These  flows  occur  in  a  surface  normal 
to  SI  and  S2,  which  is  not  considered  in  the  basic  Quasi  3D  approach. 

In  the  following,  the  different  contributions  shown  on  figure  2.1  will  be 
reviewed  separately  (sections  2.1  to  2.6).  Two  different  approaches  for  the 
determination  of  the  secondary  flows  are  reviewed  in  section  2.7.  Section 
2.8  discusses  radial  mixing  processes  and  section  2.9  is  a  brief  review  of 
experimental  work  in  this  field. 


2.1  Passage  vortex  (contributions  1,3,4) 

The  passage  vortex  is  a  well  known  secondary  flow  component  in 
turbomachine  flows.  Its  mechanism  has  been  investigated  a  first  time  by 
Carter  and  Cohen  for  uniform  flows  with  thin  boundary  layers  [26]  and  by 
Squire  and  Winter  for  nonuniform  flows  [27].  It  was  further  investigated  by 
Hawthorne  [12]  and  A.G.  Smith  [13]  for  cascades  and  by  Dixon  [14]  for 
rotating  blade  rows. 

The  passage  vortex  is  basically  a  streamwise  vorticity  component  caused  by 
the  deflection  of  a  non-uniform  mainstream  velocity.  It  can  be  obtained 
from  the  conservation  of  momentum  in  a  steady,  inviscia  incompressible  fluid 
[26,27] 

(2.11 

V  A  (W  A  w)  -  0 

The  invi3Cid  equation  (2.1)  is  often  applied  over  the  whole  passage, 
including  the  viscous  boundary  layers.  In  a  different  approach  (see  2.7) , 
the  flow  is  considered  as  a  superposition  of  (i)  a  main  inviscid  primary 
flow,  and  ( i i )  viscous  boundary  layer  flows  which  are  treated  through 


1  Main  passage  vortex 

2  Trailing  shed  and  filament  vortex 

3  End  wall  boundary  cross  flow 

4  End  wall  boundary  layer  passage  vortex 

5  Tip  clearance  flow 

6  Profile  boundary  layer  radial  flow 

7  Wake  radial  flow 

8  Radial  diffusion 


Row  non-uniformities 
Spanwise  changes  in  circulation 
Boundary  layer  overturning 
Boundary  layer  overturning 
Tip  clearance  pressure  difference 
Centrifugation 
Centrifugation 
Turbulence 


Figure  2.1 


Secondary  flows 


integral  boundary  layer  methods.  In  the  latter  case,  a  distinction  is  to  be 
made  between  a  'main  passage  vortex'  caused  by  an  overall  non-uniform  inlet 
flow  as  predicted  by  a  through-flow  computation,  and  an  'end-wall  boundary 
layer  passage  vortex'. 


2.2  Trailing  edge  vorticity  (contribution  2) 

The  vortex  sheet  leaving  the  trailing  edge  of  the  blade  has  two 
components  :  the  trailing  filament  and  the  trailing  shed  vorticities.  The 
shed  vorticity  is  due  to  the  changes  in  circulation  along  the  span,  the 
filament  vorticity  is  due  to  the  stretching  of  vortex  filaments  and  is  a 
result  of  non-uniformities  in  the  inlet  flow.  Circulation  changes  can  be 
due  either  to  an  overall  non-free  vortex  behaviour  of  the  passage  flow, 
either  to  specific  changes  in  blade  loading  inside  the  end-wall  boundary 
layers. 

These  components  were  investigated  a  first  time  by  Hawthorne  [28]  and  lame 
and  Marsh  [29]  on  the  basis  of  inviscid  flow  models. 
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2.3  Tip  clearance  flows  (contribution  5) 

Tip  clearance  flows  contribute  to  the  secondary  flow  motion  in  an 
important  way.  They  affect  the  passage  vorte-  by  modifying  the  flow  turning 
and  the  trailing  edge  vortex  through  a  change  in  blade  loading.  Smith  [11] 
proposed  an  inviscid  approach  where  tip  clearance  and  leakage  secondary 
flows  are  taken  into  account.  Several  empirical  constants  were  introduced. 
In  an  approach  where  boundary  layers  flows  are  considered  separately,  tip 
clearance  effects  can  be  included  in  the  integral  boundary  layer  equations 
[3]  (see  2.7). 


2.4  Profile  boundary  layer  radial  flows  (contribution  6) 

Boundary  layers  along  turbomachine  blades  are  three  dimensional  and 
contribute  to  the  radial  convection.  Both  experimental  and  theoretical 
information  about  these  radial  flows  is  poor.  Adkins  and  Smith  [8] 
estimated  these  radial  flows  by  solving  the  local  conservation  of  radial 
momentum  along  the  blade,  introducing  two  empirical  constants.  The  radial 
velocity  was  found  to  be  a  direct  consequence  of  the  boundary  layer  velocity 
gradient.  Their  development  is  inviscid,  in  the  same  way  as  was  done  for 
the  passage  vortex.  Whitfield  and  Keith  [ 30 ]  used  the  same  base  but 
included  the  blade  loading  effects  by  means  of  the  loss  correlations  of  Koch 
and  Smith. 

Viscous  3D  boundary  layer  applications  on  turbomachine  blades  are  almost 
inexistant.  More  experimental  information  has  recently  been  made  available 
by  Laksminara> ana  and  Popovski  [24]. 


2.5  Wake  radial  flows  (contr ibution  7) 

The  wake  can  be  considered  as  an  important  element  in  the  total  radial 
flow  analysis.  Peak  values  in  radial  velocity  are  observed  in  the  near 
wake,  where  very  low  velocities  may  induce  a  high  centrifugation.  Adkins 
and  Smith  [8]  and  Whitfield  and  Keith  [30]  treated  the  wake  by  writing  the 
mainstream  and  radial  momentum  conservation  laws  along  the  wake  center.  The 
total  amount  of  radial  flow  is  found  by  modeling  the  wake  profiles.  Since 
symmetric  wake  models  are  used,  there  is  no  ability  to  simulate  asymmetric 
radial  wake  profiles.  Hah  and  Lakshminarayana  [31]  made  full  turbulent  3D 
computations  in  asymmetric  wakes  with  radial  flows.  These  computations  are 
in  good  agreement  with  detailed  experiments,  but  their  cost  may  be  in 
disproportion  to  the  gain  in  useful  information. 


2.6  Radial  diffusion  (contribution  8) 
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above  radial  flows  contribute  to  radial  mixing  by  a  convection 
Gallimore  and  Cumpsty  [6,7]  showed  that  at  least  in  their  test 
the  radial  energy  redistribution  is  due  mainly  to  turbulent 
They  based  their  conclusion  on  the  facts  (i)  that  the  observed 
isotropic  [6]  and  (ii)  that  the  measured  turbulent  energy 
is  almost  fully  compensated  by  turbulent  energy  diffusion,  which 
convection  phenomena  are  of  less  importance  [7]. 


2.7  Two  approaches  for  secondary  flows 

From  the  the  above  discussion  on  secondary  flows,  it  follows  that  a  part 
of  them  originate  in  three  viscous  layers  : 

(i)  the  end-wall  boundary  layers 

(ii)  the  profile  boundary  layers 

(iii)  the  wakes 

The  cross  flows  in  these  layers  may  be  estimated  in  two  ways  : 

(i)  solving  local  inviscid  equations  over  the  whole  span,  including  the 
viscous  layers 

(ii)  exclude  the  viscous  layers  from  the  inviscid  solution  and 
solving  them  separately  by  means  of  integral  boundary  layer 
and  wake  equations 

The  first  approach  has  been  referred  to  in  sections  2.1  to  2.4  [8,11,12, 

1 3, 1 ^  .26,27,28,29, 30, 31  ] .  A  large  amount  of  such  inviscid  secondary  flow 
approaches  is  reviewed  in  the  AGARD  conference  proceedings  214  (1977)  about 
'Secondary  Flows  in  Turbomachines '  [25]. 

In  the  integral  boundary  layer  and  wake  approach,  secondary  flows  are 
predicted  in  a  less  detailed  way,  whereas  viscous  effects  can  be  taken  into 
account.  Secondary  flows  are  characterised  by  some  integral  cross  flow 
thicknesses  which  quantify  a  total  amount  of  secondary  mass  and  momentum 
flow. 

Early  3D  end-wall  boundary  layer  (EWBL)  equations  were  proposed  by  Cooke  and 
Hall  [15]  and  Horlock  [16].  Smith  [17]  showed  the  presence  of  variations  in 
blade  loading  inside  the  boundary  layers,  which  can  affect  the  secondary 
flow  in  an  important  way,  and  he  showed  the  importance  of  the  tip  clearance 
effects.  Mellor  &  Wood  [9]  and  Horlock  &  Perkins  [IS]  took  these  effects 
into  account  through  the  introduction  of  force  defects.  De  Ruyck  &  Hirsch 
[3]  extended  this  theory  by  introducing  (i)  new  correlations  for  the  force 
defects  and  (ii)  profile  model  equations  suited  for  turbomachine  EWBL  flows, 
which  allow  the  reconstruction  of  complete  spanwise  pitch  averaged  cross 
flow  profiles.  Only  two  empirical  constants  which  are  extensively 
calibrated  are  present  in  this  method. 

Detailed  3D  profile  boundary  layer  (PBL)  and  wake  equations  are  less  profuse 
and  as  a  f  ,'rst  step  in  the  present  project,  3D  boundary  layer  and  wake 
equations  were  developed  in  fully  curvilinear  coordinates.  This  development 
is  applied  in  a  coherent  way  for  EWBL  as  well  as  PBL  and  wakes.  Recent 
theoretical  developments  of  full  3D  boundary  layer  equations  are  summarised 
in  [25].  This  refence  discusses  the  problem  of  separations,  a  problem  which 
is  not  yet  treated  in  the  present  report,  but  which  is  to  be  included  in  the 
future  since  it  appears  from  the  available  test  cases  that  the  trailing  edge 
separation  plays  an  important  role  in  the  near  wake,  including  the  induction 
of  radial  flows. 

Although  inviscid  theories  are  important  for  the  understanding  of 
secondary  flows,  these  theories  are  practically  not  applied  in  multistage 
design  codes,  except  for  radial  flows  along  profiles  and  in  wakes.  The 
present  proposal  therefore  seeks  for  a  complete  integral  approach  for  the 
secondary  flow  description.  Integral  secondary  flow  results  are  next  used 
as  input  for  the  description  of  the  radial  mixing  process. 


2,8  Radial  mixing  process  [8] 

Adkins  &  Smith  [8]  proposed  a  diffusion  type  mixing  equation  which  can  be 
applied  for  the  mixing  of  any  flow  parameter.  The  'diffusion'  or  'mixing' 
coefficient  is  determined  from  the  knowledge  of  the  secondary  flow,  which  in 
their  case  is  obtained  from  Smith's  inviscid  secondary  flow  approach.  Some 
applications  in  [8]  show  the  importance  and  the  correct  behaviour  of  this 
mixing  theory.  More  recently,  Gallimore  and  Cumpsty  [7]  presented  much 
simpler  computations  based  on  the  turbulent  mixing  only.  Whitfield  and 
Keith  [30]  made  a  more  clear  distinction  between  convection  and  diffusion 
mechanisms.  The  did  not  consider  the  convection  as  being  isotropic,  as  was 
done  by  Adkins  through  the  use  of  a  diffusion  type  equation.  They  rather 
traced  the  flow  particles  and  obtained  non-isotropic  flow  property 
migrations  in  the  spanwise  direction. 


The  three  above  approaches  are  validaded  by  experiments,  although  the 
mechanisms  they  use  are  very  different.  This  validation  is  however  obtained 
through  the  calibration  of  empirical  constants,  which  probably  compensates 
the  shortcomings  of  these  approaches.  In  order  to  clarify  this  situation, 
one  of  the  scopes  of  the  present  work  will  be  to  determine  the  relative 
importance  of  mixing  through  convection  and  mixing  through  turbulent 
diffusion.  The  diffusion  type  equation  of  Adkins  &  Smith  will  be  applied  in 
the  present  approach  and  any  flow  property  <j>  will  be  radially  redistr ibuted 
through 


dtp  3  2ip 

3m  €  3 n'J 


(2.1  ) 


where,  combining  convection  and  diffusion, 
wn 

e  =  — j—  d0  +  turbulent  diffusion  coefficient  (2.2) 

w 

pitch  m 

w  and  w  denote  the  meridional  and  radial  velocities,  9  the  angular 
m  n 

coordinate.  Equation  (2.1)  does  not  include  mixing  losses,  it  only 
expresses  the  mixing  process.  It  can  easily  be  solved  through  a  numerical 
scheme.  The  energy  conservation  is  obtained  through  a  Neumann-type  boundary 
condition  which  expresses  that  there  is  no  enthalpy  or  entropy  flux  through 
the  end-walls.  Mixing  losses  can  eventually  be  added  through  source  terms 
in  equation  2.1. 


2 . 9  Experimental  work 

Besides  theoretical  developments,  experimental  data  can  be  found.  A 
classical  test  case  is  the  Larguier  Rotor  [5]  tested  at  ONERA.  Test  data 
are  available  downstream  of  a  rotor  blade  row,  showing  a  radial  flow  profile 
which  can  be  interpreted  as  a  superposi t ion  of  the  different  secondary  flow 
contributions  [8].  Dring  et  al  [*0  made  extensive  experiments  at  four 
stations  downstream  of  a  rotor  blade  row,  at  several  radial  positions,  and 
at  four  different  flow  rates.  Other  wake  data  are  available  from  the 
Pennsylvania  State  University  [32,33]-  Two  three-stage  compressors, 
differing  primarily  in  aspect  ratio,  were  designed  and  tested  at  Pratt  ^ 
Whitney  Aircraft  under  NASA  sponsorship  [3*1.35].  The  effects  of  mixing  ar*3 
seen  to  be  very  strong  for  both  configurations  (see  figure  1.1 '. 
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Less  detailed  data  are  available  for  the  profile  boundary  layer  radial 
flows.  Recently,  data  were  made  available  by  Lakhsminarayana  and  Popovski 
[24].  Extensive  data  are  presented  at  both  pressure  and  suction  side  of  a 
compressor  rotor  blade,  at  five  radial  positions  and  three  axial  stations. 


3-  approach 

In  the  following  sections  it  will  be  discussed  in  which  way  the  different 
radial  flows  (figure  2.1)  will  be  investigated.  There  is  an  essential 
difference  with  the  Adkins  and  Smith  [8]  and  Whitfield  and  Keith  [30  ] 
approaches  through  the  use  of  integral  methods,  reducing  the  amount  of 
empirical  constants.  The  approach  which  will  be  used  for  the  radial  mixing 
process  will  be  based  on  equations  (2.1)  and  (2.2).  The  present  report  is 
limited  to  the  radial  flow  prediction. 

The  strategy  which  is  used  is  illustrated  on  figure  3.1 •  Four  steps  can 
be  considered  in  the  global  flow  computation  : 

(i)  Through  flow  computation 

(ii)  Computation  of  end-wall  boundary  layers 

(iii)  Computation  of  the  passage  radial  flows 

(iv)  Computation  of  profile  boundary  layers  and  wakes 

All  the  equations  are  expressed  in  a  meridional  coordinate  system  m,n,u 
which  is  shown  on  figure  3-2.  The  secondary  flows  are  determined  in  a 
surface  normal  to  the  SI  and  S2  surfaces,  that  is  surfaces  of  constant 
meridional  coordinate  m. 

3.1  Through  flow  computation 

The  computation  of  the  radial  effects  needs  input  from  the  main  flow 
computations,  or  eventually  from  experiment.  In  practice  a  Quasi  32 
computation  is  performed  which  yields  the  spanwise  passage  averaged  profiles 
of  all  the  flow  properties.  Amongst  these  are  the  meridional  ana  angular 
velocity  components  w^  and  w  ,  and  the  unmixed  temperature  profiles.  The 
velocity  profiles  are  input  for  the  next  step  in  the  calculation  procedure. 
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figure  3.2  :  meridional  coordinate  system 

3-2  End-wall  boundary  layer  computation 

The  EWBL  calculation  method  a3  presented  in  [3]  can  directly  be  used.  In 
this  method  passage  averaged  EWBL  equations  are  fully  developed  in  the 
curvilinear  meridional  coordinate  system,  including  all  curvature,  Coriolis 
and  centrifugal  effects.  Appropriate  EWBL  velocity  profiles,  developed  for 
the  particular  case  of  a  turbomachine,  are  used. 

The  inclusion  of  the  boundary  layer  equations  yield  passage  averaged 

velocity  profiles  w  and  w  where 

m  u 


Vrf  *  w  -  w 
m  m  m.ewbl 


W  a  W  *  W 

u  u  u.ewbl 


(3.1  ) 


(3.2) 


Wm  ewbl  and  wu  ewbl  are  the  veiocity  defects  present  in  the  boundary  layers. 
According  to  [3],  the  resulting  passage  averaged  secondary  flow  in  the 
surface  m»cte  is  given  by 


w  ,  »  w  -  w  tgS 

u. sec. ewbl  u  m 


where  tgS 
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This  secondary  velocity  is  zero  in  case  of  a  collateral  EWBL  flow  (8  = 
implicitly  takes  into  account  the  effects  of  : 

-  viscosity 

-  tip  clearance 

-  passage  and  trailing  edge  vorticities  caused  by  EWBL  velocity  defects 
The  EWBL  calculation  procedure  is  summarized  in  section  4. 


3.3  Passage  radial  flow  computation 


The  radial  flows  inside  the  passage  are  found  from  the  incompressible 
inviscid  vorticity  equation 


(W.V).oj  -  (co.^).W  =  0 


(3.4) 


Since  the  secondary  flow  is  searched  in  surfaces  m=cte,  only  the  meridional 
component  of  the  secondary  vorticity  is  relevant 


.  3rw  ,  3w 
1  u  1  n  . 

/.)  —  i  —  —  -  a  (J 

m  r  3n  r  30 


(3.5) 


From  the  previous  steps,  some  information  is  available  about  the  passage 
average  of  the  velocity  w^,  whereas  w^  is  the  searched  velocity  component. 

The  idea  of  the  present  approach  is  to  passage  average  equations  (3.4) 
and  (3*5).  This  yields  following  equations  (section  5) 


m  _  r  r  i-  m 

—  *  F  1  u  ,  w  ,  w  ,  - — 

n  m’  1  nJ+  u’  3n 


(3.6) 


where 
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(3.7) 


The  passage  averaged  velocity  w^  can  be  decomposed  as 


w  "  w  +  w  +  vv 

u  u  u.sec.ewbl  u.sec.inv 


(3.8) 


[w  ]+  is  de  difference  in  radial  velocity  between  pressure  and  suction  side 
of  the  passage  (see  figure  3-').  Equation  (3-8)  superposes  the  velocities 
from  the  through  flow  computation  (section  3-1).  the  EWBL  flow  (section  3.2) 


bmw  iwnw  iwira  wrm vwavrwvivxvwvmvvwvm  vn  w 1 j*  ?  vrao  '.»tj.  r.»  ■■■■»  ?  m  ».w  «>  ’yv^Tvyyv 


and  main  passage  secondary  flow. 

Equations  (3.6)  to  (3.8)  contain  four  unknowns  :  u>  ,  [w  ]  ,  w  and 
-  m  n  +  u 

w  .  .  The  continuity  in  the  secondary  flow  surface  is  used  as  fourth 

u • sec • 1 n v 

equation  to  obtain  closure.  Some  pitchwise  modeling  of  the  secondary  flow 
velocities  is  required  at  this  point,  since  pitchwise  information  is  lost 
through  the  passage  averaging.  A  linear  behaviour  for  w^  and  a  quadratic 
behaviour  for  w^  yields  [9] 


9  r  8 

r—  I  W  I  =  —  W 

9n  L  n  +  s  u.sec 


(3.9) 


From  equations  (3.6)  to  (3.9),  it  is  possible  to  reconstruct  the  radial 
velocities  including 

-  all  EWBL  effects  through  w 

u.sec-ewbl 

-  main  passage  effects  through  a^w^ 

-  main  circulation  changes  through  3nr'wu 

The  only  assumptions  made  consist  in  some  pitchwise  modelling  of  the 
secondary  flow  velocities,  which  results  in  an  approximation  in  equation 
(3.9). 


3.4  Blade  profile  boundary  layers  and  wakes 

These  flows  originate  mainly  from  the  centrifugation  and  coriolis  effects 
in  blade  profiles  and  blade  wakes.  The  corresponding  3D  boundary  layer 
equations  automatically  include  these  effects.  In  the  present  wake 
approach,  pressure  and  suction  side  of  the  wake  are  treated  separately  and 
the  concept  of  'peak  velocities'  is  not  used.  The  reason  for  this  is  that 
wake  profile  may  be  very  asymmetric.  Pressure  and  suction  side  radial  flows 
may  even  have  opposite  signs  (Dring  data  [4]).  An  asymmetric  3D  wake  model 
equation  is  introduced. 

When  this  last  step  is  performed,  radial  flows  are  reconstructed  anywhere 
in  the  passage  flow  and  the  radial  mixing  process  through  radial  convection 
can  be  started. 
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4.  END-WALL  BOUNDARY  LAYER  COMPUTATIONS 


The  used  method  originates  from  Mellor  and  Wood's  approach  [9]  where 
rigorous  integral  boundary  layer  equations  were  written  for  the  end-wall 
boundary  layers,  including  force  defect  thicknesses.  In  this  work, 
simplified  assumptions  such  as  constant  shape  factor  or  skin  friction  were 
made  and  the  equations  were  integrated  from  inlet  to  outlet  of  a  blade  row. 
This  approach  was  extended  by  De  Ruyck  &  Hirsch  [3]  to  include  velocity 
profile  models  and  correlations  for  secondary  flow  and  tip  clearance  effects 
through  the  so  called  defect  forces. 

The  End-wall  boundary  theory  presently  used  has  been  presented 
extensively  in  [2]  and  is  briefly  summarized  below. 

The  basic  equations  are  written  in  the  meridional  coordinate  system  m,n,u 
(fig  3.2).  'c'  denotes  an  absolute  velocity  and  'w'  a  relative  velocity 
component . 

The  boundary  layer  momentum  and  entrainment  equations  are 
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Eqs  (4.1)  to  (4.3)  are  written  in  the  absolute  reference  system, 
boundary  layer  thicknesses  in  these  equations  are  defined  as 
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The  entrainment  rate  E  is  given  by 
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The  mainstream  and  cross  flow  shear  stresses  are  obtained  from 


Ts  "  Cf  2' 


t.  =*  t  tge 
t  s  w 
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The  'kinematic'  thicknesses  6*  and  8  .  are  defined  as 
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Eqs  (4.4)  to  (4.19)  are  written  in  the  wall  reference  system.  The 
thicknesses  0,0,6*  and  6*  and  the  parameters  b  and  n  have  the  same 


Iw*’’ 


values  in  both  absolute  and  relative  systems  whereas  the  wall  skewing  angles 


and  e  are  related  by  the  jump  relations 
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The  blade  mainstream  defect  force  is  found  from 
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The  tranverse  lift  defect  is  found  from 
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in  absence  of  a  tip  clearance  and  from 
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in  presence  of  a  tip  clearance.  The  constants  k  and  k'  are  given  by 


k  *  3  • 


k'«  .5 


(4.25) 


The  cross  flow  thickness  is  defined  through 
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Eqs  (4.21)  to  (4.26)  are  written  in  the  blade  reference  system. 

The  basic  eqs  (4.1),  (4.2)  and  (4.3)  are  integrated  in  the  meridional 
direction  using  a  fifth  order  Runge-Kutta  method  with  as  complementary 
relations  eqs  (4.4)  to  (4.26).  The  amount  of  secondary  flow  is  given  by  the 
pitch-averaged  cross  momentum  and  displacement  thicknesses,  which  will 
used  as  input  in  the  next  radial  flow  computations  (section  5). 


5.  MAIN  PASSAGE  INVISCID  SECONDARY  FLOWS 


The  ideas  of  the  present  section  have  been  given  in  section  3-3.  Passage 
averaged  secondary  flow  momentum  and  continuity  equations  will  oe  set  up,  in 
such  a  way  that  radial  velocities  will  be  obtained  including  the  EW3L 
effects,  main  passage  effects  and  main  circulation  changes. 

The  quantities  which  will  be  determined  are  (figure  5.1) 

(i)  the  inviscid  passage  averaged  cross  flow  component  w^  sgc 
(ii)  the  difference  in  radial  velocity  between  pressure  and  suction  side  [w 

[w  ]+  contains  a  contribution  from  the  end-wall  boundary  layers.  It  does 
not  contain  the  blade  profile  boundary  layer  and  wake  radial  flows  which 
will  be  treated  in  sections  6  and  8. 


Figure  5.1  :  Passage  radial  flow  computation 


5. 1  Secondary  flow  momentum  equation 

The  incompressible  inviscid  vorticity  equation  can  be  written  as 

*  ■*  ■*  ■*  -*  -*  (5.1) 

(W.V).U)  -  U.V).W  =  0 

The  vorticity  component  in  the  n  direction  (see  figure  3-2)  will  directly  be 
assumed  to  be  zero  which  means  that  the  blade  to  blade  flow  is  considered  to 
be  a  potential  flow. 


integration  of  eq  (5.2)  from  pressure  to  suction  side  yields 

-  [  w  ]  =  -  t-  rw  "  l  w  ]  +  (5.3) 

s  L  mJ  +  r  9m  u  s  L  uJ  + 

[w  ]+  represents  the  difference  in  meridional  velocity  between  suction  and 
pressure  side  of  the  blade  (fig  5.1).  s  denotes  the  pitch.  Taking  the 
following  boundary  condition  into  account 


eq  (5.3)  becomes 
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The  projection  of  equation  (5.1)  in  the  direction  u  is  of  no  concern. 
The  projection  in  direction  m  yields  (see  appendix  A) 
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R„  are  the  radii  of  curvature  of  the  m,n,u  lines  (see  appendix  A). 
Expressing  that 
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one  finds 
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Pitch  averaging  of  this  equation  yields 
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Second  order  fluctuations  have  been  neglected.  n  is  the  radial  blade 
inclination  angle  wi 
equ  (5.8)  reduces  to 


inclination  angle  which  can  be  assumed  to  be  very  small.  Assuming  w^  <<  w^ 
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The  second  term  in  the  r.h.s.  arises  from  the  presence  of  a  mainflow 

velocity  gradient  (<ou)  in  presence  of  a  flow  deflection  ([wm]+,  eq  (5.9)). 
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This  last  contribution  results  from  changes  in  circulation  along  the  blades, 
since  from  (5.3)  and  (5.9) 
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Equation  5.10  is  rigorous  when  the  radial  velocity  is  antisymmetric  with  u, 
as  will  be  assumed  below  (see  third  step  on  figure  3.2). 

The  averages  of  to  and  to  are  given  by  (metric  coefficients  are  constant 
mu 

with  x,.) 
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5.2  Superposition  of  secondary  flows 


Passage  averaged  quantities  have  been  denoted  with  a  ~  in  eqs  (5.3)  to 
(5.13).  Since  the  EWBL  flows  are  computed  separately,  care  must  be  taken 
with  the  interpretation  of  each  of  these  quantities.  They  are  to  be 
interpreted  as  follows  : 

The  meridional  velocity  is  equal  to  value  detained  from  the  through  flow 
computation. 


w  =  w 
m  m 


(5. 1  **) 


It  does  not  include  the  EWBL  velocity  defect  since  it  has  been  decided  to 
compute  the  viscous  layers  separately. 

The  pitchwise  velocity  contains  an  ’inviscid’  cross  flow  component 

w  .  which  arises  from  the  main  non-uniformities  of  the  flow  and  which 

u.sec. inv 

is  still  to  be  determined.  At  this  point,  it  is  decided  to  superpose  the 
cross-flow  contribution  from  the  EWBL  (eq  (4.7))  on  w^.  In  this  way,  fhe 
secondary  flow  information  from  the  EWBL  computation  is  considered  as  input 
for  the  present  computation,  in  the  same  way  as  for  w  .  The  rec irculat ions 
arising  from  the  EWBL  which  extend  over  the  whole  passage  area  are  hence 
automatically  included. 
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(5.15) 


5 . 3  Secondary  flow  continuity  equation 


Three  unknown  variables  appear  in  equations  (5.9)  to  (5.15) 


:  the  passage  averaged  meridional  secondary  vcrticity 


the  variation  in  radial  velocity  over  the  passage 
the  cross  flow  arising  from  main  flow  non-uniformities 


[w  ] 
Z  n 


u.sec.mv 


Closure  can  be  obtained 
continuity  equation,  written 
A) 


by  means  of  the  meridional  incompressible 
for  the  secondary  flow  velocity  W  c  (appendix 
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It  has  been  assumed  that  the  secondary  flow  can  be  defined  in  the  surface  of 
constant  meridional  coordinate  [33  and  hence  w  =0.  Integration  of  eq 
(5.16)  over  one  half  pitch  yields 
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(5.17) 


where  w  (m)  is  the  mid-pitch  cross  flow.  Equation  (5.17)  cannot  be  used 
u.sec 

without  some  spanwise  remodeling  of  the  secondary  flow  components. 
According  to  Mellor  and  Wood  [9],  the  radial  flow  will  be  assumed  to  be 
linear  and  the  pitchwise  flow  to  be  quadratic  with  the  passage  coordinate  u, 
reducing  eq  (5.17)  to 
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5.1*  Practical  solution  procedure 

The  secondary  flow  can  be  solved  from  equations  (5.9),  (5.12),  (5.13), 
(5.15)  and  (5.18).  In  practice,  equations  (5.12),  (5.15)  and  (5.1 3 )  are 
combined  to  (neglecting  curvature  for  simplicity) 
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(5.19) 


Eq  (5.19)  is  a  one  dimensional  second  order  equation  with  as  boundary 
conditions  a  zero  radial  flow  at  hub  and  casing.  It  can  easily  be  solved 
with  a  central  difference  scheme  and  inversion  of  a  tridiagonal  matrix. 


Application  in  a  straight  cascade  (r>>)  would  reduce  eq  (5. 1 9 )  to 
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which  can  easily  be  solved  analytically  when  approaching  the  r.h.s.  by  a 
constant.  Following  analytical  solution  can  be  found  in  the  case  of  a 
constant  r.h.s. 
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where 


a  =  2/2/s 
h  =  half  span 
n  =  0  at  half  span 


(5.22) 


Eq  (5.21)  is  a  simple  result.  Is  has  been  applied  to  the  Larguier  Rotor 
test  case  and  yields  a  good  agreement  close  to  that  obtained  by  solving  the 
full  equation  (5.19)  (figure  9.10). 


6.  BLADE  PROFILE  BOUNDARY  LAYERS 


The  blade  profile  calculation  is  based  on  an  adaptation  of  the  EWBL 
calculation  [3].  A  2D  adaptation  has  been  introduced  in  the  code  Q3DFL0 
described  in  [ 1  ] ,  showing  a  correct  behaviour  of  the  basic  models  used.  In 
the  present  section  the  3D  development  of  the  profile  boundary  layer  is 
presented,  which  finally  will  yield  the  radial  flow  along  the  blade 
profiles.  The  theory  is  developed  in  the  meridional  coordinate  system  which 
is  shown  on  figure  3.2  and  discussed  in  appendix  A. 

The  radial  quantities  will  are  determined  from  the  the  present  section  are 
indicated  on  figure  6.1. 


Figure  6.1  :  radial  flows  from  profile  boundary  layers 


6. 1  Momentum  equations  in  curvilinear  coordinates 


In  the  following,  the  flow  variables  denoted  with  a  carat  ~  are  those 
obtained  from  the  quasi  3D  computations.  The  integral  boundary  layer 
approach  seeks  for  the  differences  between  these  variables  and  the  real  flow 
variables. 

The  momentum  equations  will  be  derived  taking  wall  curvature  and  variable 
inviscid  flow  parameters  into  account. 

Both  real  and  inviscid  flow  satisfy  momentum  conservation.  Subtraction 


of  both  conservation  laws  yields 


p(W7)W  -  p(  W$)W  ♦  2toA  C  pW-pi? )  -  -  $(p-p)  *  $(-t) 


(6.6) 


Through  the  use  of  inviscid  flow  parameters  as  reference  ones,  radial 


equilibrium  is  taken  into  account  inside  the  P3L. 


The  m  and  n  projections  of  the  integral  boundary  layer  momentum  equation 
are  found  to  be,  after  considerable  manipulation  (appendix  B)  : 
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where  following  thickness  definitions  are  found 
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The  ±  in  equations  (6.7)  and  (6.8)  indicate  right  or  left  hand  coordinates, 
The  +  sign  is  selected  when  taking  w^  always  positive  in  the  rotation  sense, 


The  thicknesses  which  appear  in  eqs  (6.7)  and  (6.8)  originate  from 


variations  in  velocity  w.  These  terms  are  of  second  order  since  these 
variations  are  in  general  small  and  their  contribution  will  not  be 


considered.  Neglecting  X  _  reduces  these  equations  to,  after  some 


manipulations 
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All  the  variables  are  written  in  the  right  handed  blade  reference  coordinate 
system  and 


U  =  wr  in  a  rotor 
U  ■  0  in  a  stator 


For  all  approximations  made  when  obtaining  eqs  (6.13)  and  (6. 1*0  it  is 
refered  to  appendix  B. 


Eqs  (6.13)  and  (6.14)  are  to  be  written  in  the  blade  reference  system 
since  they  are  constructed  from  the  steady  Navier-Stokes  equations. 
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6.2  Entrainment  and  skin  friction 


The  entrainment  equation  is  constructed  from  the  continuity  as  follows 
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Integration  of  eq  (6.15)  with  respect  to  dx^  yields 
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Inverting  integral  and  derivative,  with  3h  /3x  =  0  and  h  dx  =  rd9 
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Division  by  h  h  and  rearranging  leads  to 
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which  is  easily  rearranged  to 
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The  last  term  is  related  to  the  entrainment  rate  E  as  follows,  s  being  the 
streamline  direction 
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Eq  (6.20)  is  written  in  the  relative  coordinate  system.  The  entrainment 
rate  is  given  as  an  empirical  function  of  Head's  shape  factor  through  : 

E  -  0.O306/(H*-3.)-653  <6'2,) 


The  validity  of  this  empirical  relation  is  questionable  in  presence  of 
curved  walls  but  no  corrections  to  this  equation  are  found  in  the  open 
literature.  Since  it  is  however  shown  that  curvature  has  a  strong  effect  on 
the  turbulence  properties  of  the  flow,  it  is  believed  that  a  correction 
should  be  introduced  in  eq  (6.21).  No  attempts  are  made  in  the  present  work 
to  correct  eq  (6.21)  and  this  should  be  considered  as  a  limitation  in 
presence  of  strongly  curved  walls. 


The  estimation  of  the  mainstream  skin  friction  occurs  through  the  use  of 
a  non  dimensional  skin  friction  coefficient  generally  denoted  as  and 
defined  through 
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Ludwig-Tillman ' s  relation  is  generally  accepted.  This  relation  expresses 
as  function  of  the  Reynolds  number  and  shape  factor  H  and  is  originally 
given  by 


C_  -  0.246  Re  exp(-1.56  H) 
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As  for  the  entrainment  rate,  Sumner  and  Shanebrook  [19]  analysed 
compressibility  effects  on  skin  friction  and  showed  the  validity  of  the 
following  equation 
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S3 

where  Eckert's  reference  temperature  T*  and  the  kinematic  shape  factor  H  are 
introduced.  The  reference  temperature  T*  for  adiabatic  walls  is  given  by 

T*  Y-1  2 

—  =  1  *  0.72  r  -z-  M  (6.26) 
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e 

where  r  is  the  recovery  factor. 

Experiments  show  that  the  effect  of  curvature  on  the  skin  friction  is  not  to 
be  neglected.  These  experiments  however  indicate  that  the  Ludwig-Tillman 
equation  eq  (6.25)  is  still  valid.  For  different  types  of  curved  boundary 
layer  flows  the  increase  in  friction  and  in  shape  factor  due  to  curvature 
are  related  by 


(6.27) 


This  behaviour  is  correctly  predicted  by  the  Ludwig-Tillman  correlation , 
since  from  eq  (6.25)  it  follows  that 


1  .56 


(6.28) 


Eqs  (6.27)  and  (6.28)  are  practically  equivalent  for  current  values  of  H. 


6.4  Profile  models 


Velocity  profiles  are  introduced  in  order  to  close  the  system  of 
equations.  Velocity  profiles  and  boundary  layer  thicknesses  have  been 
defined  in  meridional  coordinates  (figure  6.2).  In  this  way,  the 
orthogonality  of  the  basic  equations  is  conserved  and  all  the  curvature, 
Coriolis  and  centrifugal  contributions  are  found  in  a  simple  way.  In 
analogy  with  the  EWBL  profile  models  (section  4),  following  model  equations 
are  introduced  : 
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where  s  lines  are  the  projection  of  streamlines  on  an  axisymmetric  stream 
surface  and  where  n  lines  are  normal  to  it. 


Streamwlse  profile 
written  in  function  of 
u  coordinate 


u 


radial  profile 
written  in  function  of 
u  coordinate 


figure  6.2  :  profile  boundary  layer 
velocity  profiles 


The  parameters  b  and  n  are  correlated  through  the  following  relation  [2,3, 
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The  introduction  of  the  model  equations  in  the  different  thicknesses  is 

detailed  in  appendix  C. 


6.5  Solution  procedure 

The  following  differential  equations  are  considered  : 

(i)  the  meridional  momentum  equation  (6.13) 

(ii)  the  radial  momentum  equation  (6.1*0 

(iii)  the  entrainment  equation  (6.19) 

Extra  relations  are 

(i)  the  entrainment  rate  equation  (6.21) 

(ii)  the  friction  law  equation  (6.23) 

(iii)  the  profile  model  equations  (6.29)  and  (6.30), 
applied  in  the  definitions  (6.9)  to  (6.12) 


Radial  momentum  and  displacement  thicknesses  are  obtained,  from  which  radial 
profiles  can  be  reconstructed  through  the  radial  velocity  model  equation 


7.  WAKE  ANALYSIS 


The  way  used  by  both  Adkins  and  Smith  [8]  and  Whitfield  and  Keith  [30]  to 
analyse  the  wake  is  to  consider  the  momentum  equation  at  the  wake  centerline 
and  to  assume  a  symmetric  wake  profile.  Radial  flows  are  characterised  by 
the  'peak'  radial  velocity  at  the  wake  center.  In  the  present  wake 
approach,  pressure  and  suction  side  of  the  wake  are  treated  separately  and 
the  concept  of  ’peak  velocities'  is  not  used.  The  reason  for  this  is  that 
wake  profiles  may  be  very  asymmetric.  Pressure  and  suction  side  radial 
flows  may  even  have  opposite  signs,  as  can  be  seen  from  the  experiments  of 
Dring  [4], 

The  contribution  to  the  radial  flows  which  is  obtained  from  the  present 
section  is  shown  on  figure  7.1. 


Figure  7.1  :  wake  radial  flows 


7 . 1  Wake  profiles 


The  definition  of  the  wake  profiles  occurs  in  the  same  way  as  has  been 
done  in  section  6.^  for  the  blade  profiles.  The  basic  model  equations  are 
(see  figure  7.2) 
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The  function  f  is  shown  to  be  close  to  a  probability  function  given  by 


which  can  be  approached  by  [22] 

f(n)  =  1  -  6n2  +  Sn3  -  3n“  (7.*0 

The  wake  is  in  general  considered  to  be  symmetric,  which  in  reality  is  not 
the  case.  For  this  reason,  pressure  and  suction  side  of  the  wake  will  be 
considered  separately  and  a  correction  in  its  central  part  will  be 
introduced  as  follows,  e.g.  for  the  pressure  side  (  +  ) 
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where  the  correction  function  g  chosen  as 
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The  +  or  -  signs  indicate  pressure  or  suction  sides.  The  overbars  indicate 
mean  values  such  as 


b  =  b~  (7.8) 
These  profiles  are  illustrated  on  figure  7.2. 

The  function  g  corrects  the  profile  in  the  neighbourhood  of  the  centerwake 
and  makes  a  bridge  between  pressure  and  suction  sides.  A,t  n=0,  the 
velocities  reduce  to  (for  both  sides) 
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Figure  7.2  :  wake  profile  model 


For  small  values  of  n0.  pressure  and  suction  side  of  the  wake  conserve  their 


identity,  at  higher  values  both  sides  are  mixing  together  and  the  wake  shape 
becomes  more  symmetric.  High  second  derivatives  are  present  around  the 


centerline  and  the  asymmetry  gives  rise  to  a  diffuse  mixing  between  pressure 


and  suction  side  of  the  wake.  Hence  n0  can  be  approximated  by 
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where  Am  is  the  distance  from  the  trailing  edge. 


7.2  Integral  wake  momentum  equations 


Since  pressure  and  suction  side  of  the  wake  are  considered  separately, 


two  sets  of  equations  are  used.  These  equations  can  be  shown  to  be  exactly 
the  same  as  the  blade  momentum  equations  (6.13)  and  (6.14).  Due  to 
asymmetry,  the  friction  terms  in  eqs  (6.13)  and  (6.14)  do  not  disappear,  as 
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long  as  b  *  b  and  tge  *  tge  .  The  friction  is  given  by  (with  n=0) 
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The  pressure  and  suction  side  frictions  are  related  by 
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As  a  consequence,  the  friction  tends  to  decrease  the  difference  between 
momentum  thicknesses  at  pressure  and  suction  side. 


From  equations  (7.4)  to  (7.11)  it  is  possible  to  determine  all  the  integral 
thicknesses,  defined  through  eqs  (6.9)  and  (6.11).  The  corresponding 
expressions  for  the  wake  thicknesses  and  their  derivatives  are  listed  in 
appendix  D. 


7.3  Local  wake  decay  equations 


The  decay  of  the  center  wake  velocity  is  determined  from  the  center  wake 
momentum  equation  which  can  be  written  as  (see  section  6.2) 
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(7.15) 
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Assuming  small  and  ~  0  at  the  wake  center,  the  meridional  component 

can  be  found  from  (see  appendix  B) 
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Using  eqs  (7.16)  to  (7.18),  equation  (7.15)  reduces  to 
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Applying  eq  (7.5),  and  neglecting  the  curvature  term  for  a  while,  one  can 
find,  for  n=0  at  pressure  side 
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The  term  in  g"(0)  gives  a  discontinuity  in  the  second  derivative  when  going 
from  pressure  to  suction  side.  Taking  the  average  over  both  sides  and 
equation  (7.4)  for  f  (f"(0)=-l2)  one  can  write,  for  b 
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where  a  mean  value  for  5  is  defined  as 
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Expressing  the  momentum  equation  (7.15)  in  n=0,  only  one  equation  is 
obtained  for  both  sides.  This  equation  describes  the  mean  wake  velocity 


defect . 


In  order  to  obtain  more  information  about  pressure  and  suction  si^e 
mixing,  a  relation  is  searched  between  b  and  b  .  This  relation  is  obtained 
by  writing  eq  (7.15)  in  nx,  where  nx  is  such  that 

1  =  fx»  g,  ;7-22' 


n,  is  a  point  close  to  the  wake  center  but  which  is  not  influenced  by  the 
other  wake  side.  In  this  way,  eq  (7.15)  reduces  to,  for  the  pressure  side 
and  suction  side 
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and,  taking  the  difference,  with  Ab  »  b  -  b 
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and  taking  eqs  (7.^)  and  (7.7)  for  f  and  g  yields,  with  eq  (7.22) 
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In  a  very  similar  way,  by  considering  the  projection  of  eq  (7.15)  in  the  n 
direction,  one  can  find 
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Eqs  (7.26)  and  (7.27)  close  the  system  of  wake  equations,  the  only  unknown 
left  is  the  eddy  viscosity  coefficient  v  which  is  obtained  from 


(7.28) 


v  =  K5bw 

m 

The  coefficient  K  is  to  be  calibrated. 

7 . 4  Practical  solution  procedure 

Six  independent  variables  are  present  in  the  wake  equations  : 

5,6  :  pressure  and  suction  side  wake  thickness 

+  **• 

b  ,  b  :  pressure  and  suction  side  wake  velocity  defect 

+ 

e  ,  e  :  pressure  and  suction  side  skewing  angle 

Closure  is  obtained  with  the  following  differential  equations  : 

pressure  and  suction  side  eqs  (6.13)  : 

integral  meridional  momentum  expressed  at 
pressure  and  suction  side 

sum  of  pressure  and  suction  side  eq  (6.14)  : 

sum  of  the  integral  radial  momentum  expressed  at 
pressure  and  suction  side 

equation  (7.21)  :  local  meridional  momentum  equation  along  wake  center 

equation  (7.26)  :  difference  between  local  meridional  momentum 

equations  expressed  at  pressure  and  suction  side 

equation  (7.27)  :  difference  between  local  radial  momentum 

equations  expressed  at  pressure  and  suction  side 

The  two  radial  momentum  equations  (6.14)  could  be  considered  at  pressure  and 
suction  side  separately,  making  equation  (7.27)  redundant.  In  practice, 
numerical  instability  was  however  observed  when  using  the  radial  momentum 
equations  separately  at  pressure  and  suction  side.  This  instability  is  due 
to  singularities  in  the  Jacobian  matrix  which  relates  the  momentum 
thicknesses  with  the  wake  parameters  b,  6  and  tge.  in  particular,  it  was 
found  that  the  profile  model  equations  are  not  very  sensitive  to  the 
difference  Atge  and  the  radial  momentum  equation  (6.14;  has  therefore  been 
expressed  over  the  whole  wake.  This  equation  mainly  controls  the  mean 


skewing  tge  and  dissappears  from  the  system  of  equations.  Equation 

(7.13)  for  x  is  not  used, 
n 

New  parameters  which  appear  in  these  equations  are  the  integral  wake 
momentum  and  displacement  thicknesses,  the  central  wake  shear  stress  ts  ,  a 
mixing  length  n0  and  a  diffusion  coefficient  v.  These  extra  parameters  are 

found  from 


equations  (7.5)  and  (7.6)  :  wake  velocity  models  which  allow  the 

determination  of  the  integral  wake  thicknesses 

equation  (7.12)  :  stress  from  the  model  equations 

equation  (7.11)  :  model  for  the  mixing  length  v0 

equations  (7.28)  :  model  for  the  diffusion  coefficient  n 

¥  -  ♦  - 

The  meridional  profiles  are  modeled  by  5  ,  6  ,  b  ,  b  .  These  variables  are 

the  main  parameters  of  the  two  meridional  integral  equations  (6.13)  and 

meridional  local  equations  (7.21)  and  (7.26).  The  radial  profiles  are 

¥  - 

modeled  by  the  skewing  angles  e  and  e  .  The  difference  Atge  i3  controlled 
by  the  extra  equation  (7.27). 
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8.  TRAILING  EDGE  AND  NEAR  WAKE 


The  wake  equations  as  discussed  in  section  7  are  not  directly  applicable 
in  a  region  very  close  to  the  trailing  edge,  mainly  because  of  the  presence 
of  strong  pressure  gradients  which  are  not  considered  in  the  momentum 
equations  (6.13)  and  (6.14).  A  detailed  description  inside  this  region  is 
not  considered,  since  the  complexity  of  the  computation  is  not  in  proportion 
to  the  gain  in  information.  Direct  jump  relations  between  the  momumtum  and 
displacement  thicknesses  at  trailing  edge  and  downstream  wake  are  therefore 
searched. 

These  relations  are  based  on  work  of  Kool  and  Hirsch  [20].  The  pressure 
recovery  from  trailing  edge  to  the  downstream  wake  center  is  proposed  as 
[20] 


C  =  0.3  (1  -  56  C  ) 
P  Q 


(8.1  ) 
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Application  of  these  relations  in  the  wake  momentum  equations  yields  [20] 
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It  i3  clear  that  this  relation  cannot  be  valid  for  the  complete  range  of  the 

ratio  r  /e^  .  For  very  small  trailing  edge  radii,  one  should  in  the  limit 
te  te 
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The  wake  equations  as  discussed  in  section  7  are  not  directly  applicable 
in  a  region  very  close  to  the  trailing  edge,  mainly  because  of  the  presence 
of  strong  pressure  gradients  which  are  not  considered  in  the  momentum 
equations  (6.13)  and  (6.1*0.  A  detailed  description  inside  this  region  is 
not  considered,  since  the  complexity  of  the  computation  is  not  in  proportion 
to  the  gain  in  information.  Direct  jump  relations  between  the  momumtum  and 
displacement  thicknesses  at  trailing  edge  and  downstream  wake  are  therefore 
searched. 

These  relations  are  based  on  work  of  Kool  and  Hirsch  [20].  The  pressure 
recovery  from  trailing  edge  to  the  downstream  wake  center  is  proposed  as 
[20] 

C  =  0.3  (1  -  56  C  )  (8.1 ) 
p  q 

where 


,v->j 


Experiments  [20]  indicate  that  eq  (8.5)  is  valid  from  rte/6te  *  10  on.  For 
lower  ratio's  a  simple  linear  interpolation  between  eqs  (8.5)  and  (8.6)  is 
proposed  through 

9  =  0t  ( 1  +0.0i)5r,  )  +  .01  r .  /8.  (8.7) 

R  te  te  te  te 

The  momentum  thickness  is  not  sufficient  as  information  to  determine  the 
wake  profile  downstream  of  the  near  wake  region.  According  to  Kool  and 
Hirsch  [20],  one  should  take  b  close  to  1  according  to  the  existance  of  a 
stagnation  point  inside  the  near  wake,  b  should  however  not  be  regarded  as 
the  real  velocity  defect,  but  rather  as  a  parameter  in  the  model  equation  eq 
(7.1).  This  parameter  probably  deviates  from  its  physical  meaning  in  the 
near  wake  and  it  controls  the  whole  outer  wake,  where  the  profile  boundary 
layer  history  is  still  strongly  present.  For  these  reasons,  it  is  proposed 
to  keep  the  physical  thickness  6  as  a  constant  value  in  the  jump  over  the 
trailing  edge,  rather  than  imposing  an  arbitrary  value  for  b. 

A  third  jump  relation  is  required  for  the  radial  flow  thicknesses. 
Assuming  no  particularities  in  the  radial  pressure  gradient  in  the  near 
wake,  one  can  assume  the  radial  momentum  to  be  constant. 

In  conclusion,  the  jump  relations  can  be  written  as 

9  -  f ( 0  .  ,r  ) 

mm  mm,te  te 
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te 

9*0. 
nn  nn,te 

where  eq  (8.5)  or  (8.7)  apply  for  the  first  relation 


9.  RESULTS 


Three  test  cases  were  considered  at  present  :  an  isolated  airfoil  wake 
which  is  two  dimensional  [23],  the  Larguier  rotor  [5]  and  the  Dring  rotor 
[4]. 


9.1  Isolated  airfoil  [23] 

The  isolated  airfoil  case  is  used  to  test  the  wake  equations  in  the 
two-dimensional  case  for  both  symmetric  and  asymmetric  flow  behaviour.  The 
airfoil  is  a  NACA-0012  airfoil,  experiments  are  available  at  0°,  5°  and  10° 
angle  of  attack.  Results  are  shown  on  figures  9.1  to  9.6.  Figures  9.1,  9.3 
and  9.5  compare  the  experimental  and  the  computed  wake  displacement  and 
momentum  thicknesses,  on  figures  9.2,  9.4  and  9.6  the  corresponding  wake 
velocity  profiles  are  shown.  The  wake  thicknesses  are  compared  separately 
for  the  pressure  side,  suction  side  and  total  wake.  The  asymmetric  results 
show  the  correct  transfer  of  momentum  from  the  thicker  suction  side  to  the 
thinner  pressure  side  of  the  wake  (figures  9.3  and  9.5).  The  transfer  of 
momentum  occurs  only  through  the  shear  stress  is  which  is  computed  from 
equation  (7.12),  indicating  a  correct  estimation  of  n0  through  equation 
(7.11).  The  coefficient  K  is  found  to  be  0.01. 

The  wake  profiles  (figures  9.2,  9.4,  9.6)  show  the  correct  behaviour  of 
the  model  equations  (7.4)  to  (7.7)  for  downstream  positions  of  5,  10  and  20 
%  chord  distance  from  the  trailing  edge.  At  shorter  distances  0.7  and  2  % 
chord,  a  better  agreement  cannot  be  found  since  the  wake  profiles  are  still 
closer  to  the  boundary  layer  profile  shapes  than  to  the  Gaussian  wake  shape. 
Attempts  to  define  a  transition  region  in  the  near  wake  were  undertaken,  but 
the  gain  in  information  is  not  in  proportion  to  the  increase  of  complexity 
required.  The  transition  region  between  profile  boundary  layer  and  wake 
computations  is  rather  to  be  bypassed  by  the  use  the  jump  relations  8.8. 


9.2  ONERA  rotor  [5] 

The  ONERA  rotor  has  been  considered  as  a  test  case  for  radial  mixing 
computations  by  several  authors  [8,30].  Results  are  presented  near  the  mean 
diameter  at  30  %  axial  chord  distance  downstream  of  the  trailing  edge.  The 
blades  have  radially  constant  solidity  and  no  twist,  giving  rise  to  a  a 
substantial  central  secondary  flow. 
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The  results  presented  in  this  section  are  obtained  by  a  full  computation, 
without  any  input  from  the  experimental  data.  The  QNERA  rotor  was  treated 
by  a  Quasi  3D  code  which  yields  all  the  input  required  to  run  the  EWBL,  main 
passage  secondary  flow.  Profile  Boundary  Layers  and  wake  computations. 

The  new  elements  which  are  tested  at  this  level  are  the  profile  boundary 
layer  computations  which  result  in  the  entry  conditions  of  tr.e  waxe 
computations,  and  main  secondary  flow  equations  5.9  and  5.18. 

The  result  obtained  at  trailing  edge  is  shown  on  figure  9.7.  It  is  to  be 
observed  that  the  radial  velocities  observed  in  the  Larguier  rotor  are 
mainly  due  to  its  non-free  vortex  design.  The  contribution  from  EWBL  and 
profile  boundary  layer  secondary  flows  is  very  small,  as  can  be  seen  from 
figure  9.8  where  the  vorticity  was  taken  zero. 

Figures  9.9  and  9.10  compare  the  experimental  wake  profiles  with  computed 
results.  The  differences  between  computation  and  experiment  lie  within  the 
differences  observed  in  [5]  when  different  types  of  anemometers  are  used. 

9.3  Dring  rotor  [4] 

Extensive  measurements  were  made  by  Dring  et  al  [4]  downstream  of  a 
single  rotor  blade  row.  Velocity  measurements  are  available  at  four 
downstream  stations,  at  several  radii.  Measurements  were  made  at  four 
different  flow  rates. 

It  is  the  aim  of  the  present  work  to  predict  these  results  through  a  total 
computation  (Quasi3D,  EWBL,  PBL,  central  secondary  flow,  wake).  The  profile 
boundary  layer  computations  however  predict  trailing  momentum  thicknesses 
which  are  far  too  small  (7  to  10  times  the  observed  values)  and  at  present 
this  problem  has  not  yet  been  solved.  The  underestimation  is  due  either  to 
an  inaccurate  estimation  of  the  velocity  gradients  along  the  profiles, 
either  to  the  presence  of  separations  near  the  trailing  edge,  which  may 
cause  a  momentum  loss  which  is  much  larger  than  the  one  predicted  through  eq 
(8.5). 

At  present,  the  testing  of  the  Dring  rotor  is  therefore  limited  to  the  3D 
wake  computation,  using  the  first  experimental  station  data  as  input. 


Results  are  presented  at  two  different  flow  rates  ($  »  0.65  and  j>  =■  0.85) 
on  figures  9.11  and  9.12.  These  figures  compare  the  obtained  meridional  and 


radial  flow  velocities  with  Dring's  data,  at  three  radii.  The  lowest  flow 
rate  is  close  to  the  stall  limit,  whereas  the  higher  one  gives  a  certain 
amount  of  secondary  flow  due  to  non-free  vortex  behaviour.  A  quite  good 
agreement  is  in  general  observed.  The  value  for  the  dispersion  coefficient 
K  in  eq  (6.23)  was  however  to  be  increased  considerably  for  toe  test  case 
clcse  to  the  stall  limit,  which  can  be  explained  only  by  an  important 
increase  in  turbulence. 

A  strong  radial  flow  is  observed  experimentally  inside  the  wake  close  to  the 
stall  limit.  Such  a  strong  radial  flow  is  not  found  when  applying  the 
profile  boundary  layer  equations.  It  is  probably  due  to  separations  in  the 
neighbourhood  of  the  trailing  edge,  causing  a  stronger  wake  centrifugation 
at  this  place,  as  stated  by  Dring  et  al. 

The  total  radial  recirculation  due  to  the  non-free  vortex  behaviour  at  the 
higher  flow  rate  appears  to  be  more  important  then  the  wake  radial  flow 
itself,  indicating  that  strong  radial  mixing  effects  are  to  expected  in 
off-design  conditions. 


CONCLUSIONS 


The  main  tools  required  for  the  prediction  of  the  radial  flow 
distributions  in  a  compressor  stage  'nave  been  developed  in  the  present 
report.  This  theory  is  based  on  the  superposition  of  the  different 
contributions  from  end-wall  boundary  layers,  profile  boundary  layers,  wakes 
and  main  passage  secondary  flows.  All  the  viscous  aspects  are  treated  by  an 
integral  approach  where  velocity  profile  models  are  introduced.  The 
inviscid  secondary  flow  behaviour  is  obtained  from  pitch  averaged  vorticity 
and  continuity  equations. 

The  behaviour  of  the  asymmetric  wake  and  of  the  central  inviscid  secondary 
flows  are  well  verified  by  experiments.  The  behaviour  of  the  EWBL  flows  are 
treated  elsewhere  [2,3].  The  profile  boundary  layers  is  applied  in  the 
Larguier  test  case,  where  a  full  computation  has  been  done.  No  direct 
comparison  is  made  for  the  Profile  boundary  layers  however. 

Strong  radial  flows  are  induced  by  a  non-free  vortex  behaviour, 
indicating  possible  strong  radial  mixing  effects  in  off-design  conditions. 
Strong  radial  flows  are  also  observed  in  the  wake,  close  to  the  stall  limit, 
these  strong  radial  flows  are  not  predicted  by  the  profile  boundary  layers 
equations  in  the  present  state  of  the  project. 

In  the  present  state,  the  theory  which  has  been  developed  seems  to  be 
able  to  predict  the  radial  flow  in  a  reliable  way  when  no  important 
separations  occur,  as  can  be  expected  at  design  conditions.  In  this 
situation  radial  recirculations  should  be  due  mainly  to  the  EWBL  and  PSL 
cross  flows.  Through  the  approach  which  has  been  used,  all  the  secondary 
flow  aspects  are  implicitly  or  explicitly  present  in  the  radial  flow 
prediction,  including  effects  such  as  a  tip  clearance  through  the  EWBL 
theory. 
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FUTURE  WORK 


The  main  part  of  the  radial  mixing  theory  lies  in  the  prediction  of  the 
radial  flow  intensity,  which  implies  a  total  knowledge  of  the  secondary 
flows  in  a  turbomachine. 

Each  aspect  of  this  secondary  flow  prediction  is  to  be  tested  separately, 
which  has  not  been  done  for  the  profile  boundary  layers  equations  yet. 
Recently  detailed  Profile  Boundary  Layers  experimental  data  have  been  made 
available  [24,25]  and  a  separate  testing  of  the  profile  boundary  layers  will 
be  performed. 

Once  the  radial  flows  known,  the  mixing  process  itself  can  be  described. 
The  main  difficulty  is  to  balance  the  contributions  radial  mixing  due  to 
convection  and  mixing  due  to  turbulent  diffusion  when  calibrating  the  two 
contributions.  This  will  be  analyzed  in  the  final  report. 
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APPENDIX  A  :  MERIDIONAL  COORDINATE  SYSTEM 


A. 1  Some  notions  of  curvature 

In  the  present  appendix  some  notions  of  curvature  are  reviewed  and 
applied  to  the  vector  operations  such  as  divergence  and  rotation. 

It  can  be  shown  that  in  each  orthogonal  coordinate  system,  a  set  of 
coordinates  x^x^.x^  with  unit  vectors  e^e^.e^  can  be  defined  in  such  a  way 
that  all  the  surfaces  x^cte,  x2=cte,  x^=cte  intersect  in  three  orthogonal 
curves.  x1 ,  x2  and  x^,  as  well  as  their  dimensions,  are  a  priori  unknown 
f unctions ,  except  for  special  cases  of  orthogonal  coordinates.  In  a 
cylindrical  coordinate  system  for  example,  x^.x^.x^  are  found  to  be  x,r,8 
with  as  dimensions  L,L,1. 

The  length  61  along  the  coordinate  line  6x^  is  by  definition  given  by 


61 .  »  h . 6x  . 
l  li 


(  A  .  1  ) 


where  hi»h2,h3  are  the  lame  or  metric  coefficients,  which  are  also  a  priori 
unknown  functions  (e.g.  the  cylindrical  coordinate  system  :  h1»h2=1,  h^-r) 


V 

Figure  A.1  :  Elementary  volume  in  curvilinear  coordinates 


The  gradient  operation  can  easily  be  expressed  in  a  curvilinear 
coordinate  system,  since  from 


it  easily  follows  that 


7<j> 
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When  applying  the  nabla  operator  to  a  vector  or  a  tensor,  derivatives  of  the 
unit  vectors  e^  however  appear  which  are  non-zero  due  to  curvature.  Hence 
correction  terms  in  general  appear  when  applying  equations  of  motion  in 
curvilinear  systems.  In  order  to  determine  these  corrections,  the  general 
definition  of  the  nabla  operator  will  be  introduced  [Milne]  (fig  A.1) 
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7.X  =  lim  -  S  n.X  dS 
,,  V  s 
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(A.  4) 


where  X  denotes  as  well  a  scalar  as  a  vector  or  a  tensor.  Choosing  an 
elementary  volume 
easily  found  that 


elementary  volume  dV  with  dimensions  h.dx.,  and  applying  eq  A. 4,  it  is 
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From  eq  A. 5  it  fellows  that  e.g.  for  e. 
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Hence,  eqs  A. 5  and  A. 6  can  be  rewritten  as  (sum  over  i  and  j) 
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where  t .  =  r .  .e  . 
1  ij  J 


In  eqs  (A. 8)  and  (A. 9),  the  correction  terms  due  to  curvature  appear  as 
functions  of  derivatives  of  the  unit  vectors. 


Since  the  length  corresponding  to  5x^  is  h  fix..,  it  can  easily  be  seen 
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from  fig  A. 2  that  for  example 


1  3el 


h29x2 


1  - 

"  R  S2 
n21 


(A. 10) 


where  R  is  the  radius  of  curvature  of  the  plane  curve  obtained  by 
intersection  of  the  surface  x^cte  with  the  plane  (e^Se^,  or  (e1  .e^)  (fig 
A. 2).  According  to  the  theorem  of  Meusnier,  this  radius  of  curvature  can  be 
related  to  R2,  the  principal  radius  of  curvature  of  x2  through  (fig  A. 2) 
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Figure  A. 2  :  Radii  of  curvature  in  curvilinear  coordinates 


From  eqs  A. 7  and  A. 10  repeated  for  each  R  it  can  be  shown  that  e.g. 


for  e, 
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Hence,  if  all  the  radii  R  are  known  functions,  the  correction  terms  in  eqs 
A. 8  and  A. 9  are  known,  while  the  functions  and  h.  can  -principally-  be 
determined  from 
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A. 2  The  meridional  coordinate  system 

The  meridional  coordinate  system  is  illustrated  on  fig  3.1.  In  general 

curvature  is  present  in  both  circumferential  and  meridional  directions.  The 

corresponding  coordinates  are  m,n,u  where  m  denotes  the  meridional 

direction,  n  the  direction  normal  to  the  end-wall  and  u  is  the  pitchwise 

direction,  m,  n  and  u  are  defined  as  lengths.  The  independent  coordinates 

are  denoted  as  x  ,  x  and  x  and  the  corresponding  metric  coefficients  as 

m  n  u 

h  ,  h  ,  h  . 

m  n  u 

The  direction  u  is  taken  positive  in  the  direction  normal  to  the  profile 
walls,  while  n  is  always  taken  positive  from  the  hub  wall.  Hence,  if  the 
pressure  side  coordinate  system  is  right  handed,  the  suction  side  coordinate 
system  will  be  left  handed,  and  vice-versa. 

According  to  the  theorem  of  Dupin,  the  system  of  coordinates  m,n,u  is  an 
orthogonal  coordinate  system,  since  the  m  and  u  lines  are  lines  of  principal 
curvature.  A  rigorous  development  in  non-orthogonal  streamline  coordinates 
is  not  considered. 

The  different  radii  of  curvature  of  the  coordinate  lines  along  the  blade 
profiles  are  found  to  be  (fig  3.1) 
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According  to  section  A.1,  the  knowledge  of  these  radii  allows  the 
application  of  any  vector  operation  in  the  meridional  coordinate  system. 


From  eqs  (A. 14),  (A. 19)  and  (A. 20),  it  follows  that 


x  =  e 
u 


h  =  r 
u 


(A. 21  ) 


(A. 22) 


h  ,  h  ,  x  and  x  cannot  be  determined  explicitly  for  an  arbitrary  wall 
m  n  m  n 

shape.  Due  to  axisymmetry,  it  however  follows  from  eqs  (A. 14),  (A. 16)  and 
(A. 18)  that 


(A. 23) 


It  is  to  be  observed  that  eqs  (A. 15)  to  (A. 20)  are  valid  in  both  pressure 
and  suction  side  m,n,u  coordinate  systems  which  are  respectively  right  and 
left  handed  or  vice-versa.  Care  must  be  taken  when  chosing  the  sign 
convention  of  o. 
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APPENDIX  B  :  MOMENTUM  EQUATIONS  IN  CURVILINEAR  COORDINATES 


In  the  present  appendix,  the  momentum  equations  are  reconstructed  taxing 
wall  curvature  and  variable  i.nviscid  flew  parameters  into  account.  Ir,  the 
following,  the  flow  variables  denoted  with  a  carat  *  are  those  obtained  from 
the  quasi  3D  computations.  The  integral  boundary  layer  approach  seeks  for 
the  differences  between  these  variables  and  the  real  flow  variables. 

Both  flows  must  satisfy  continuity,  subtraction  of  both  conditions  yields 
(in  the  blade  reference  system) 


$.(pW  -  pW)  =  0 


(B.  1  ) 


According  to  eq  (A. 5),  eq  (B.1)  can  be  written  as 
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=  0 


Eq  (B.2)  will  be  integrated  over  the  boundary  layer.  The  integration  of  eq 
(B.2)  yields 
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where 


w  =  0  at  x  =  0 
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W  =  W  at  x  =  x  r 
u  u  6 


Using  eqs  (A.1),  (A. 21)  and  (A. 22),  eq  ( B . 3 )  can  be  rewritten  as 


(3.5) 
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The  subscript  w  stands  for  'wall'.  It  is  to  be  observed  that  the 
derivatives  have  been  put  outside  of  the  integral  sign  in  eq  (3.5).  This 
may  be  done  since  the  function  to  be  integrated  becomes  zero  at  x  (5). 


The  momentum  equation  3.3  for  variable  inviscid  flows  is  to  be  rewritten 


p(W$)W  -  p(W$)W  +  2uA(pW-pW) 


p(p-p)  +  p(t-t) 


In  eq  B.6,  the  difference  in  momentum  between  the  inviscid  and  the  real  flow 
is  equated.  Through  the  use  of  inviscid  flow  parameters  as  reference  ones, 
radial  equilibrium  is  taken  into  account  inside  the  EWBL.  In  the  following 
the  four  contributions  in  eq  B.6  will  be  analysed  separately. 


B.1  Inertial  contribution 


The  mainstream  projection  of  (WV)W  is  found  to  be 
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Hence,  taking  h  h  h  times  eq  B.7,  the  e  projection  of  the  inertial 
mnu  m 

contribution  of  eq  B.6  can  be  rewritten  as 
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1 /R  is  zero  due  to  axis/mmetry.  The  last  three  terms  in  eq  B.8  are 


correction  terms  due  to  wall  curvature. 


Integration  of  the  third  term  in  the  Ihs  of  eq  B.8  with  respect  to  x 


over  the  boundary  layer  yields 
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The  first  term  of  B.9  reduces  to 
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and,  according  to  eq  5.19  to 


-  /  w  r—  (pw  -  pw  )h  h  dx 
mw  3x  m  m  n  u  u 
m 


-  /  w  t—  (pw  -  pw  )h  h  dx 
mw  3x  n  n  m  u  u 
n 


According  to  the  continuity  equation  B.5,  the  second  term  of  3 . J  can  be 


written  as 
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With  eqs  B.9.B.11  and  B.12,  and  according  to  eqs  A.1  and  A. 14,  eq  3, 
easily  be  rearranged  to 
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In  a  similar  way,  the  spanwise  projection  of  the  inertial  term  is  found 
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B.2  Coriolis  contribution 


The  second  contribution  in  eq  B.6  will  now  be  integrated  over  the 
boundary  layer.  This  contribution  which  represents  the  coriolis  effect  can 
be  rewritten  as 
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Eq  B . 1 6  may  be  written  in  both  pressure  and  suction  side  coordinate  systems, 
the  upper  and  lower  signs  correspond  respecively  to  right  and  left  handed 
coordinate  systems.  According  to  eqs  B . 1 5  and  B.16,  the  m  and  n  projections 
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of  the  coriolis  term  are  found  to  be  respectively 
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As  for  the  first  contribution  of  eq  B.6,  eqs  B.17  and  B.18  are  multiplied  by 
hmhn^u  and  next  integratecl  with  respect  to  dx^over  the  boundary  layer  giving 
respectively 
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B.4  Shear  stress  contribution 

The  last  contribution  which  is  finally  to  be  analysed  is  the  shear  stress 
contribution.  According  to  eq  A. 6,  this  contribution  can  be  written  as  (sum 
over  j ) 
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The  meridional  projection  of  eq  B.21  is  found  to  be,  according  to  A. 10,  A. 12 
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Written  as  a  defect  and  integrated  over  the  boundary  layer  with  respect  to 

dx  the  mainstream  contribution  of  the  shear  stress  terms  is  easily  found  to 
u 


S  h  h  h  e  . V( x-x)dx  = 
m  n  u  m  u 


/  t—  ( t  -  t  )h  rd0  +  /  t—  (x  -  t  )h  rd0 
3x  mm  mm  n  3x^  nm  nm  n 

m  n 


[h  h  (t  -  i  )  ] 

L  m  n  am  um  Jo 


(B.23) 


+  r  f  "nn~  JB  -  -Jg ll  +  T.u.u_-T.^i  )  h  h  roe 

J  1  ft  R  R.  ;  m  n 


nm 


mn 


um 


only  the  central  term  is  hold,  assuming  the  other  terms  of  smaller  order  of 

magnitude.  Assuming  the  shear  stress  tangential  at  the  blade  wall,  the 

blade  shear  stress  t  can  be  written  as 

s 

x  =  1  cosg  x  +  1  sing  x  (B.2M) 
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On  the  other  hand,  x  can  be  transformed  from  the  other  shear  stresses 
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APPENDIX  C  :  INTEGRAL  BOUNDARY  LAYER  PARAMETERS 


In  this  appendix  integral  boundary  layer  parameters  are  listed  as 


functions  of  the  boundary  layer  parameters  6,  n,  b  and  tge'  .  The  thickness 

W 


6  is  not  an  independent  variable  and  it  is  shown  how  it  can  be  eliminated. 
Basic  thicknesses  are  expressed  as  direct  functions  of  these  parameters,  and 
others  are  derived  from  the  basic  thicknesses.  The  different  components  of 
the  Jacobian  matrix  are  listed.  All  the  equations  are  written  in  the 
meridional  wall  reference  system. 


The  basic  model  equations  which  are  used  are  remembered  to  be 


1  -  b  (i  -  n)n 


(c.i ) 


w  -  w  tga 

— — -  =  (i-b)tge^  (i  -  n)n 

w 

s 


(C.2) 


with  as  complementary  relations 
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The  meridional  thicknesses  are  defined  through 
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The  integrations  in  C.5  and  C.6  are  taken  along  u  lines.  Working  along 
lines  normal  to  the  blade  surface  yields  following  expressions 
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According  to  eqs  C.l  to  C.8,  it  is  easily  shown  that 
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The  other  thicknesses  are  easily  found  to  be 
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The  relation  between  the  two  sets  of  thicknesses  is  found  from 


rd0  *  6'dn/cosS 

w  »  w  cosB 
m  s 


w  »  w  sinB 
u  s 


Application  of  C.17  in  the  different  definitions  yields 
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Using  eqs  C.9  to  C .  1 8  all  the  classical  thicknesses  can  be  determined 

from  the  knowledge  of  n,  b  and  e'.  The  shape  factio  H  is  to  found  from 
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The  independent  variables  which  are  used  in  the  method  are  defined  as 
follows 


{q}  =  {  ln(-lnb)  ,  ln(n)  ,  tge'  } 

w 


( C .  20  ) 


The  corresponding  components  of  the  jacobian  matrix  J  between  the 
thicknesses  and  the  parameters  q  are  obtained  as  follows.  Differentiation 
of  eq  C.H  yields 
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where,  according  to  eqs  C.3 
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which  explains  the  choice  of  q  in  eq  C.20.  dH  is  given  by 
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The  derivatives  of  5-5*  finally  follow  from  the  definition  of  H* 
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Dif ferentiation  yields 
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The  derivatives  of  H*  are  found  to  be 
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APPENDIX  D 


INTEGRAL  WAKE  PARAMETERS 


In  this  appendix  integral  wake  parameters  are  listed  as  functions  of  the 

wake  parameters  6,  b  and  tge'.  Basic  thicknesses  are  expressed  as  direct 

w 

functions  of  these  parameters,  and  others  are  derived  from  the  basic 
thicknesses.  The  different  components  of  the  Jacobian  matrix  are  listed. 
All  the  equations  are  written  in  the  meridional  blade  reference  system. 

The  basic  model  equations  which  are  used  are 


=  1  -  b  f(n)  -  (b-b  )  g(n) 


( D .  1  ) 


h  +■  + 

•»  w  tga 

- — ^ -  =»  (1-b+)tge*f (n)  +  ((1-b)(tge  +  tga  -  tga*)  -  (1-b+)tge  )g(n)  (D.2) 


where 


f(n) 


1  - 

g(n) 


6n2  +  8n3 
=  e’n/n° 


-  3n- 


n0 


7  /vAm/w 
5  m 


(D.3) 

(D.4) 

(D.5) 


where  n0should  not  exceed  about  .5.  Thicknesses  in  streamline  coordinates 
are  defined  through 
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According  to  eqs  D.1  to  D.7,  it  is  easily  shown  that 
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using  eqs  D.3  and  D.4  for  f  and  g  yields 


/ f  dn  = 


/f2dn 


/ g  dn  -  n0 


/g2dn  - 


/fgdn  =  n0  - 


3  o  4  5 

I2n«  ♦  48n0  -  72n0 


The  meridional  thicknesses  are  found  from  (see  appendix  C) 
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Application  of  D.19  in  the  different  definitions  yields 
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